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We extract an effective strong coupling constant using low-Q 2 data and sum rules. 
^ ' Its behavior is established over the full Q 2 -range and is compared to calculations 

00 , based on lattice QCD, Schwinger-Dyson equations and a quark model. Although 

00 ■ the connection between all these quantities is not known yet, the results are sur- 

prisingly alike. Such a similitude may be related to quark-hadron duality. 

o\ . 

<^> , 1. The strong coupling constant 

*o ■ 

. A peculiar feature of strong interaction is asymptotic freedom: quark-quark 

interactions grow weaker with decreasing distances. Asymptotic freedom is 
O ' expressed in the vanishing of the QCD coupling constant, a s (Q 2 ), at large 

1^ ■ Q 2 . Conversely, the fact that a s (Q 2 ), as calculated in pQCD, becomes 

| large when Q 2 — > A.q CD is often linked to quark confinement. Since it 

_C ■ is not expected that pQCD holds at the confinement scale and since the 

condition a s (Q 2 ) — > oo when Q — > Aqcd is far from necessary to assure 
confinement 1 , it is interesting to study a s (Q 2 ) in the large distance domain. 
. Experimentally, moments of structure functions are convenient objects 

to extract a s . Among them, Y\~ n is the simplest to use. In pQCD, it is 
linked to the axial charge of the nucleon, gA, by the Bjorken sum rule: 

rr*= f\x{g{{x)-gUx)) = \9A[l-^-Z.w{^-) 2 (1) 
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where g\(g\) is the first spin structure function for the proton (neutron). 
The fit(Q 2 )/Q t 2 are higher twist corrections and become important at 
lower Q 2 . This series, usually truncated to leading twist and to 3rd or- 
der, can be used to fit experimental data and to extract a s . The higher 
twists can be computed with non-perturbative models or can be extracted 
from data, although with limited precision at the moment 3 . This imprecise 
knowledge and the break down of pQCD at low Q 2 prevent a priori the 
extraction of a s at low Q 2 . However, an effective strong coupling constants 
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was defined by Grunberg 4 in which higher twists and higher order QCD 
radiative corrections are incorporated. Eq. 1 becomes by definition: 

Tr n ^\g A [l-^.]. (2) 

D 7T 

This definition yields many advantages: the coupling constant is extractable 
at any Q 2 , is well-behaved when Q 2 — ► Aqcd, is not renormalization 
scheme (RS) dependent and is analytic when crossing quark thresholds. 
The price to pay for such benefits is that it becomes process-dependent 
(hence the subscript g\ in Eq. 2). However, as pointed out by Brodsky et 
al. 5 , effective couplings can be related to each other, at least in the pQCD 
domain, by "commensurate scale equations" . These relate, using different 
Q 2 scales, observables without RS or scale ambiguity. Thus, one effective 
coupling constant is enough to characterize the strong interaction. 

Among the possible observables available to define an effective coupling 
constant, Y\~ n has unique advantages. The generalized Gerasimov-Drell- 
Hcarn (GDH) 6 ' 7 and Bjorken sum rules predict T\~ n at low and large Q 2 , 
and T\~ n is experimentally known between these two domains. Hence, 
a Stgi can be extracted at any Q 2 . In particular, it has a well defined value 
at Q 2 =0. Furthermore, we will see that a SjSl might best be suited to be 
compared to the predictions of theories and models. 
2. Experimental determination of ct s , gi 

A measurement of T\~ n at intermediate Q 2 was reported recently 8 and was 

c 

Si' 



used to extract a.,.^ 9 . The results are shown by the triangles in Fig. 1, to 



gcther with a Sj9l extracted from SLAC data at Q =5 GeV (open square). 
Note that the elastic contribution is not included in T^~ n . 
Y\~ n is related to the generalized GDH sums: 

T p-n = _Q_fQ DH P _ GDH n^ (3) 

where a is the QED coupling constant. Hence, at Q 2 —0, T^~ n = and 

a s,gi = (4) 
At Q 2 — 0, the GDH sum rule implies: 

*r'-£i<f"> B ' -<">«•) («) 

where k p (k u ) is the proton (neutron) anomalous magnetic moment. Com- 
bining Eq. 2 and 5, we get the derivative of a s ^ gi at Q 2 =0: 

dQ 2 4g A ( M 2 M 2) - W 
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Relations 4 and 6 constrain a s , gi at low Q 2 (dashed line in Fig. 1). At large 
Q 2 , T^~ n can be estimated using Eq. 1 at leading twist and a s calculated 
with pQCD. a Si9l can be subsequently extracted (gray band). 

These data and sum rules give a SjSl (Q 2 ) at any Q 2 . A similar result is 
obtained using a model of Y\~ n and Eq. 2 (dotted line). The Burkert-Ioffc 11 
model is used because of its good match with data. 

One can compare our result to effective coupling constants extracted 
using different processes. a s ^ T was extracted from T-decay data 12 from the 
OPAL experiment (inverted triangle). It is compatible with The 
Gross-Llewellyn Smith sum rule 13 (GLS) can be used to form a s ,F 3 . The 
sum rule relates the number of valence quarks in the hadron, n v , to the 
structure function F^Q 2 ^). At leading twist, it reads: 



F?,{Q 2 , x)dx = n v 



« S (Q 2 



3.58 



a s (Q 2 



20.21 



a s (Q 2 



(7) 



We expect a Si F 3 



at high Q 2 , since the Q 2 -dependence of Eq. 1 



and 7 at leading twist are identical. The GLS sum was measured by the 
CCFR collaboration 14 and the resulting a s ,_F 3 is shown by the star symbols. 



<Z S s /7t JLab 
□ a s Jk world data 



Burkert-loffe 

GDH constraint 

I I pQCD evol. eq. 
T a, /7l OPAL 




Q (GeV) 

Figure 1. Extracted a 3>gi (Q)/tT using JLab data (up triangles), the GLS sum rule 
(stars), the world n data (open square), the Bjorken sum rule (gray band) and the 
Burkert-loffe Model. a SjT (Q) /it from OPAL is given by the reversed triangle. The 
dashed line is the GDH constrain on the derivative of a s ,g 1 /'T at Q 2 =0. 
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3. Comparison with theory 

Just like effective coupling constants extracted experimentally, there are 
also many possible theory definitions for the coupling constant and, con- 
trarily to the experimental quantities, the relations between the various 
definitions are not well known. Furthermore, the connection between the 
experimental and the theoretical quantities is not clear. Hence, the remain- 
der of this paper is to be understood as a candid comparison of quantities 
a priori defined differently, in order to see if they share common features. 

Calculations of a s using Schwinger-Dyson equations (SDE), lattice QCD 
or quark models are available. Different SDE results are shown in Fig. 2. 
The pioneering result of Cornwall 15 is shown by the blue band in the top 
left panel. The more recent SDE results from Fisher et al., Bloch et al, 
Maris and Tandy, and Bhagwat et al. are shown in top left, top right, 
bottom left and bottom left panels respectively There is a good match 
between the data and the result from Fisher et al. and a fair match with 
the curve from Bloch et al. The results from Maris- Tandy, Bhagwat et al. 
and Cornwall do not match the data. The Godfrey and Isgur curve in the 
top right panel of Fig. 2 is the coupling constant used in the framework 
of hadron spectroscopy 20 . Q 2 -behavior of coupling constants can also be 
compared regardless of their absolute magnitudes by normalizing them to tt 
at Q 2 = (These curves are not shown here). The Godfrey-Isgur, Cornwall 
and Fisher et al. <3 2 -behavior match well the data. The normalized curves 
from Maris- Tandy, Bloch et al. and Bhagwat et al. are slightly below the 
data (by typically one sigma) for Q > 0.6 GeV. 

Gluon bremsstrahlung and vertex corrections contribute to the running 
of a s . Modern SDE calculations include those 21 but it is a priori not the 
case for the a s used in the one gluon exchange term of the Godfrey and 
Isgur quark model, or for older SDE works. If so, pQCD corrections should 
be added to these calculations. The effect of those corrections (on a s , gi ) 
is given by the ratio of a SjSl extracted using Eq. 2 to a Sj9l extracted using 
Eq. 1 at leading twist. For both Eq. 1 and 2, T\~ n is given by a model 11 . 
Since model and data agree well, no strong model dependence is introduced. 
The difference between results using Eq. 1 up to A th and 5 th order is taken 
as the uncertainty due to the truncation of the pQCD series. The resulting 
a s are shown in the bottom right panel of Fig. 2. 

Finally, we can compare lattice QCD data to our results. Many lattice 
results are available and are in general consistent. We chose to compare 
with the results of Furui and Nakajima 22 , see bottom left panel in Fig. 2. 
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They match well the data. The lowest Q 2 point is afflicted by finite size 
effect and should be ignored. 

The match between our data and the various calculations might be 
surprising since these quantities are defined differently We can try to un- 
derstand this fact. Choosing " minimizes the role of resonances, in 
particular it fully cancels the A1232 contribution which usually dominates 
the moments at low Q 2 . By furthermore excluding the elastic contribution, 
we obtain a quantity for which coherent reactions (elastic and resonances) 
are suppressed and we are back to a DIS-like case in which the interpre- 
tation is straightforward. One can also possibly invoke the phenomenon 
of quark-hadron duality to explain why the extraction of a s , 9l , using a 
formalism developed for DIS 12 , seems to also work at lower Q 2 . 

4. Conclusion 

We have extracted, using JLab data at low Q 2 together with sum rules, an 
effective strong coupling constant at any Q 2 . A striking feature is its loss of 
Q 2 -dependence at low Q 2 . We compared our result to SDE and lattice QCD 
calculations and to a coupling constant used in a quark model. Despite 
the unclear relation between these various coupling constants, data and 
calculations match in most cases, especially for relative Q 2 -dependences. 
This could be linked to quark-hadron duality. 
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We extract an effective strong coupling constant using low-Q 2 data and sum rules. 
Its behavior is established over the full Q 2 -range and is compared to calculations 
based on lattice QCD, Schwinger-Dyson equations and a quark model. Although 
the connection between all these quantities is not known yet, the results are sur- 
prisingly alike. Such a similitude may be related to quark-hadron duality. 

1. The strong coupling constant 

A peculiar feature of strong interaction is asymptotic freedom: quark-quark 
interactions grow weaker with decreasing distances. Asymptotic freedom is 
expressed in the vanishing of the QCD coupling constant, a s (Q 2 ), at large 
Q 2 . Conversely, the fact that a s (Q 2 ), as calculated in pQCD, becomes 
large when Q 2 — > Aqcd ls often linked to quark confinement. Since it 
is not expected that pQCD holds at the confinement scale and since the 
condition a s (Q 2 ) — > oo when Q — > Aqcd is far from necessary to assure 
confinement 1 , it is interesting to study a s (Q 2 ) in the large distance domain. 

Experimentally, moments of structure functions are convenient objects 
to extract a s . Among them, T\~ n is the simplest to use. In pQCD, it is 
linked to the axial charge of the nucleon, qa, by the Bjorken sum rule: 

r?-" = f dx(g\(x) - gUx)) = \g A [l - - - 3.58 (^Y (1) 
Jo o 71 v 71 J 

3 4 5 00 P — 11 
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, 2 Q 

where g\ (g™) is the first spin structure function for the proton(neutron). 
The Ht(Q 2 )/Q t ~ 2 are higher twist corrections and become important at 
lower Q 2 . This series, usually truncated to leading twist and to 3rd or- 
der, can be used to fit experimental data and to extract a s . The higher 
twists can be computed with non-perturbative models or can be extracted 
from data, although with limited precision at the moment 3 . This imprecise 
knowledge and the break down of pQCD at low Q 2 prevent a priori the 
extraction of a s at low Q 2 . However, an effective strong coupling constants 
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was defined by Grunberg 4 in which higher twists and higher order QCD 
radiative corrections are incorporated. Eq. 1 becomes by definition: 

rr n = \g A [i - (2) 

This definition yields many advantages: the coupling constant is extractable 
at any Q 2 , is well-behaved when Q 2 — > Aqcd, is not renormalization 
scheme (RS) dependent and is analytic when crossing quark thresholds. 
The price to pay for such benefits is that it becomes process-dependent 
(hence the subscript gi in Eq. 2). However, as pointed out by Brodsky et 
al. 5 , effective couplings can be related to each other, at least in the pQCD 
domain, by "commensurate scale equations". These relate, using different 
Q 2 scales, observables without RS or scale ambiguity. Thus, one effective 
coupling constant is enough to characterize the strong interaction. 

Among the possible observables available to define an effective coupling 
constant, Y\~ n has unique advantages. The generalized Gerasimov-Drell- 
Hearn (GDH) 6 ' 7 and Bjorken sum rules predict T\~ n at low and large Q 2 , 
and Y\~ n is experimentally known between these two domains. Hence, 



a S}9l can be extracted at any Q . In particular, it has a well defined value 
at Q 2 =0. Furthermore, we will see that a S:9l might best be suited to be 
compared to the predictions of theories and models. 
2. Experimental determination of a s , gi 

A measurement of Y\~ n at intermediate Q 2 was reported recently 8 and was 
used to extract a S:9l 9 . The results are shown by the triangles in Fig. 1, to- 
gether with a S:9l extracted from SLAC data 10 at Q 2 =5 GeV 2 (open square). 
Note that the elastic contribution is not included in Y\~ n . 
Y\~ n is related to the generalized GDH sums: 

t p-u = _Q^_ <q DH p _ GDH n^ / 3 ) 

where a is the QED coupling constant. Hence, at Q 2 =0, Y\~ n = and 

a S}9l = TT. (4) 
At Q 2 = 0, the GDH sum rule implies: 

T '' n = S- a {GDH '- GDH ''^ z ^ { w f - J k <5) 

where k p (k„) is the proton (neutron) anomalous magnetic moment. Com- 
bining Eq. 2 and 5, we get the derivative of a St9l at Q 2 =0: 

das, gi _ 37T K 2 p 

X ^ 1 r> 1 fl " W 



dQ 2 Ag A y M 2 M 2 V 
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Relations 4 and 6 constrain a St9l at low Q 2 (dashed line in Fig. 1). At large 
Q 2 , T\~ n can be estimated using Eq. 1 at leading twist and a s calculated 
with pQCD. a St9l can be subsequently extracted (gray band). 

These data and sum rules give a Stgi (Q 2 ) at any Q 2 . A similar result is 
obtained using a model of T\~ n and Eq. 2 (dotted line). The Burkert-Ioffe 11 
model is used because of its good match with data. 

One can compare our result to effective coupling constants extracted 
using different processes. a SiT was extracted from r-decay data 12 from the 
OPAL experiment (inverted triangle). It is compatible with a. s>Bl , The 
Gross-Llewellyn Smith sum rule 13 (GLS) can be used to form a St p s . The 
sum rule relates the number of valence quarks in the hadron, n v , to the 
structure function F^(Q 2 ,x). At leading twist, it reads: 




We expect a Sj F 3 = Q? Sjgl at high Q 2 , since the <5 2 -dependence of Eq. 1 
and 7 at leading twist are identical. The GLS sum was measured by the 
CCFR collaboration 14 and the resulting a St p 3 is shown by the star symbols. 



a sgl /nJLab 
□ a s Jk world data 



Burkert-Ioffe 

GDH constraint 

I I pQCD evol. eq. 
T a, /7l OPAL 




Q (GeV) 

Figure 1. Extracted a s ,g 1 (Q)j'K using JLab data (up triangles), the GLS sum rule 
(stars), the world Y\~ n data (open square), the Bjorken sum rule (gray band) and the 
Burkert-Ioffe Model. a s , T {Q)/Tt from OPAL is given by the reversed triangle. The 
dashed line is the GDH constrain on the derivative of a Stgi /-K at Q 2 =0. 
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3. Comparison with theory 

Just like effective coupling constants extracted experimentally, there are 
also many possible theory definitions for the coupling constant and, con- 
trarily to the experimental quantities, the relations between the various 
definitions are not well known. Furthermore, the connection between the 
experimental and the theoretical quantities is not clear. Hence, the remain- 
der of this paper is to be understood as a candid comparison of quantities 
a priori defined differently, in order to see if they share common features. 

Calculations of a s using Schwinger-Dyson equations (SDE), lattice QCD 
or quark models are available. Different SDE results are shown in Fig. 2. 
The pioneering result of Cornwall 15 is shown by the blue band in the top 
left panel. The more recent SDE results from Fisher et al., Bloch et al., 
Maris and Tandy, and Bhagwat et al. are shown in top left, top right, 
bottom left and bottom left panels respectively. There is a good match 
between the data and the result from Fisher et al. and a fair match with 
the curve from Bloch et al. The results from Maris-Tandy, Bhagwat et al. 
and Cornwall do not match the data. The Godfrey and Isgur curve in the 
top right panel of Fig. 2 is the coupling constant used in the framework 
of hadron spectroscopy 20 . Q 2 -behavior of coupling constants can also be 
compared regardless of their absolute magnitudes by normalizing them to ir 
at Q 2 = (These curves are not shown here). The Godfrey-Isgur, Cornwall 
and Fisher et al. Q 2 -behavior match well the data. The normalized curves 
from Maris-Tandy, Bloch et al. and Bhagwat et al. are slightly below the 
data (by typically one sigma) for Q > 0.6 GeV. 

Gluon bremsstrahlung and vertex corrections contribute to the running 
of a s . Modern SDE calculations include those 21 but it is a priori not the 
case for the a s used in the one gluon exchange term of the Godfrey and 
Isgur quark model, or for older SDE works. If so, pQCD corrections should 
be added to these calculations. The effect of those corrections (on a St9l ) 
is given by the ratio of a St9l extracted using Eq. 2 to a s , 9l extracted using 
Eq. 1 at leading twist. For both Eq. 1 and 2, T\~ n is given by a model 11 . 
Since model and data agree well, no strong model dependence is introduced. 
The difference between results using Eq. 1 up to 4 th and 5 th order is taken 
as the uncertainty due to the truncation of the pQCD series. The resulting 
a s are shown in the bottom right panel of Fig. 2. 

Finally, we can compare lattice QCD data to our results. Many lattice 
results are available and are in general consistent. We chose to compare 
with the results of Furui and Nakajima 22 , see bottom left panel in Fig. 2. 
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They match well the data. The lowest Q 2 point is afflicted by finite size 
effect and should be ignored. 

The match between our data and the various calculations might be 
surprising since these quantities are defined differently. We can try to un- 
derstand this fact. Choosing T\~ n minimizes the role of resonances, in 
particular it fully cancels the A1232 contribution which usually dominates 
the moments at low Q 2 . By furthermore excluding the elastic contribution, 
we obtain a quantity for which coherent reactions (elastic and resonances) 
are suppressed and we are back to a DIS-like case in which the interpre- 
tation is straightforward. One can also possibly invoke the phenomenon 
of quark-hadron duality to explain why the extraction of a St9l , using a 
formalism developed for DIS 12 , seems to also work at lower Q 2 . 

4. Conclusion 

We have extracted, using JLab data at low Q 2 together with sum rules, an 
effective strong coupling constant at any Q 2 . A striking feature is its loss of 
<2 2 -dependence at low Q 2 . We compared our result to SDE and lattice QCD 
calculations and to a coupling constant used in a quark model. Despite 
the unclear relation between these various coupling constants, data and 
calculations match in most cases, especially for relative Q 2 -dependences. 
This could be linked to quark-hadron duality. 
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